Extraction of Pure Entangled States from Many Body Systems by Distant Local 
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We study the feasibility of extracting a pure entangled state of non-complementary, and potentially 
well separated, regions of a quantum many-body system. It is shown that this can indeed be 
accomplished in non-equilibrium scenarios as well as the ground state of the considered spin chain 
models when one locally measures observables such as magnetization in separated blocks of spins. 
A general procedure is presented, which can search for the optimal way to extract a pure entangled 
state through local projections. Our results indicate a connection of the projective extraction of 
entanglement to good quantum numbers of the underlying Hamiltonian. 



I. INTRODUCTION 

Much attention as been paid to entanglement, the 
purely quantum part of correlations, since the emergence 
of the quantum information science. In this context, 
quantum entanglement can be regarded as a resource 
for quantum information processors [3]. As example of 
this assertion, entanglement is believed to be the ma- 
jor resource for the acceleration capabilities of quantum 
processors in performing computation and communica- 
tion tasks such as teleportation or superdense coding. 
With this in mind, it is important to note that a typical 
quantum processor is also a quantum many body system. 
Thus, the attention that quantum information commu- 
nity has paid in the last years to quantum many body 
systems intensively studied in condensed matter physics 
has lead to consider that condensed matter phenom- 
ena can be explored from an information theoretic point 
of view, particularly studying entanglement in quantum 
many body systems [2| and using these systems for fruit- 
ful operations in quantum communications [H and com- 
putation protocols [J]. 

In considering the use of many body systems for quan- 
tum information applications, one would ideally like to 
establish pure entangled states between distant regions, 
because all quantum protocols such as teleportation p|, 
work best for pure states. There is also the additional 
technical advantage that one can use the von Neumann 
entropy of entanglement @ , which is a very good measure 
for pure states to compute the entanglement. Indeed, for 
ground states which are pure, von Neumann entropy has 
been extensively used to calculate entanglement between 
complementary parts 0], [H, @- In this context the 
amount of entanglement obtainable in useful form from 
these complementary parts for a single copy of a many- 
body state has also been discussed [10(. But for generic 
distant parts of a QMB system the state is a mixture, 
so a conversion to a pure state would be ideal. It would 
be desirable to achieve this by only local operations on 
those parts. Of course, in principle, one could use several 
copies and employ purification protocols, which require 
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FIG. 1: Localizable Entanglement (Top), i) Representation of 
a quantum many body system, ii) To localize entanglement 
between the two end parties, we perform local measurements 
on the other parties, iii) LE is the maximum amount of en- 
tanglement that on average can be localized on these end 
parties by this procedure. Projective Extraction of Pure En- 
tanglement (Bottom), i) Representation of a quantum many 
body system, ii) We want to establish a pure entangled state 
between regions A and B by means of two coarse grained pro- 
jective measurements on these regions, iii) The procedure can 
establish pure entangled states between regions A and B with 
some probability of success. 



local actions as measurements and unitary operations as 
well as classical communication between the parties, to 
achieve a smaller number of almost pure states. How- 
ever, it is highly desirable that these local operations be 
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simple (such as projections) acting on a single copy of an 
entangled state. 

In [ill ] was introduced a new measure of bipartite en- 
tanglement called Localizable Entanglement (LE) that is 
operationally defined as the maximum amount of entan- 
glement that can be localized, on average, between two 
spins of the chain by performing local measurements on 
other parties. Studies on the amount of LE as a function 
of the distance of the spins to be entangled and the influ- 
ence of anisotropic couplings were conducted in [13] .By 
definition, the LE procedure requires a fine grained ac- 
cess to all the subsystems of the many body system. This 
requirement seems rather difficult for typical situations 
where the subsystems are very close to each other. 

In this paper we present a procedure that, in a proba- 
bilistic sense, allows us to establish pure bipartite entan- 
gled states between distant regions of many body system 
by means of local projections only on these regions (Figl). 



II. PROJECTIVE EXTRACTION OF PURE 
ENTANGLED STATES 

We consider the state \<j>) of a system S defined on a 
lattice and the associated density operator po — \4>)((j>\- 
We can define projective operators Pa, Pb acting locally 
on separate regions A and B of S, where each of them 
may comprise several lattice sites. The state after local 
projection PaPb is 



= (PaPb)po(PaPb) m 
9 Tr(P A P B po) 1 ' 

By an appropriate choice of the projectors Pa, Pb 
the state pab = T r (S-{AuB))(p) Wli l be pure, and 
the probability of obtaining this state when perform- 
ing the associated selective measurement is Prob(pAB) — 
Tr(PAPBPo)- Once this pure state is established, the en- 
tanglement between regions A and B can be quantified 
by the von Neumann entropy S(pa) = — Tr(pA log 2 Pa) 
with pa = Ttb(pab)- It is convenient to introduce 

£pp{Pa, Pb) = S(p A ) Prob(p AB ) (2) 

where the subscript letters PP indicate the features 
projective and pure. This quantity both captures the 
probabilistic aspects of this procedure and quantifies the 
amount of quantum correlations in the resulting state. 
The procedure is meaningful only for a choice of the pro- 
jectors Pa and Pb that are (i) suitable for extraction 
of pure states and (ii) of rank higher than one because 
otherwise one would always obtain an unentangled state. 
We shall make these requirements more concise in Section 

ITZ3 



III. SOME SIMPLE EXAMPLES OF £ PP 

In this section we illustrate how the projective proce- 
dure can be applied in two simple examples. 

A. Supersinglets 

Supersinglets are defined as [lU 

\S ( n) = "4= € su-,s N \si,S2,--.,S N ) , (3) 

describing a highly entangled state of a number iV qu- 
dits with individual number d of degrees of freedom. The 
summation runs over all permutations si...jv of ./V-tuplcs 
of si, each of which assuming integer values [1, d] and la- 
beling the local basis states {|sz)} on site I, and e sll ... iSJV 
denotes the generalized Levi-Civita tensor. It was re- 
cently found that supersinglet states also arise as ground 
states for permutation Hamiltonians of systems where 
d = N (generalization of two spin-i coupled by a Heisen- 
berg interaction) [l4[ ■ These states are particularly well 
suited to illustrate the concept of pure entangled state 
extraction by local projections, as introduced in the pre- 
vious Section. To this end we choose a simple system 
of three qutrits (N = d = 3), (labelled A, B and C) the 
supersinglet state reads 

|*) = -^(|i) c |a> 2 , 3 > AB + |2) c |$ 3 ,iU + m c \*i,2) AB ) 

(4) 

where we defined = — By choosing, 

e.g. 

P A/B = |2)(2| + |3)(3| (5) 

after projection the qutrits A and B will be left in the 
pure state |$23) for which S($23) = 1- In a correspond- 
ing selective measurement, this will happen with proba- 
bility i. Then we find £pp = \- To hightlight require- 
ment (ii), we stress that by choosing rank-1 projector, 
e.g. 

P A = \2){2\ 

P s = |3><3] (6) 

the state after projection will be pure but separable, 
namely |2)^|3) B . Rank-3 projectors for a qutrit act triv- 
ially as the identity, and leave the subsystems A and B 
in a maximally mixed state. For general d = N one finds 

£ pp = Tr ■ 

B. Epp in a non-equilibrium scenario 

In this Section, we consider a ring of N spins interact- 
ing through the exchange Hamiltonian 

N 
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with periodic boundary conditions a% 
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or a 



x,y, 



ri rz 



where a x and a y are the Pauli a matrices. The initial 
state of the system is not an eigenstate of ((7]) and consists 
of two spin flips located at sites r\ and r 2 of the ring 



^=(|n)|r 2 ) 



(8) 



This physical situation allows temporal evolution for 
i]) giving the state for certain time t, \<f>t) 
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J2 e(s 2 - Sl )A SuS2 (t)\ Sl )\s 2 ) (9) 



where A SyS2 (t) = - f£ ri (t)f£ r2 (t)) is 

the amplitude that at time t, one of the flips being at 
site si and the other one being at site s 2 and 
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is the amplitude for one of the flips to be initially at 
r and at time t be located at s, with fj — At, d = s — r, 
f = 0.5 and e(i — j) is the sign function (e(x) — +1, x > 
0; e(x) = —1, x < 0; e(x) = 0, x = 0) and J u (x) is a Bcsscl 
function of the first kind. 

The objective of our procedure is to create a pure en- 
tangled state between two separate regions A and B of 
the ring by projecting ([9]) with two local coarse grained 
measurements which ascertain whether there is some flip 
present in certain region without revealing any informa- 
tion about its position within the region (see illustration 
in Fig. [J). This procedure discards trials in which no 
flips are found in any region or only one of the flips is 
found in one the measurement regions. The projected 
state after measurement is the result of selecting the cases 
in which we find one flip in A and the other one in re- 
gion B. This setting of our projection procedure can be 
made more concise by explicitely writting the projectors 
as P A = EjeAliXil and P b = E fce sl fc )( fc l- So, tak- 
ing into account © and (TIT))) the state after projection 
I^*)ab = PaPb\^h) can be written as 
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e(k-j)A jk (t)\j)\k) (11) 



with A jk (t) = (f n (t) N f k2 (t) N -f j2 (t) N f kl (t) N ) and also 
noting that Vj S A,k G B e(k - j) = +1. 

In order to compute the Epp resulting for our partic- 
ular choice of projectors Pa and Pb we have to compute 
the von Neumann entropy of the reduced density matrix 
of one of the regions, that in the basis of lattice sites 
(real space basis) it is easy to show that has the form, 
(pa(%/ = (j'\PA(t)\j) = EkeB^3k(t)A*, k (t) where j 

and j £ A. The probability for obtaining the state (jXTJ) 
in our procedure also easily reads as 



Prob(\V t ) AB ) 



2£ 



jeA.keB 



E 



(12) 



llllltltlllll 



II 1 III 1 1 1 III f II 



nliiiiiinii 



1 1 1 ii 1 1 1 ii 1 1 ir i 




FIG. 2: Schematic diagram of the Epp in the XY ring with 
two spin flips. The initial state of consists in two localized flips 
in an initially non-entangled state (Top). After some time 
evolution, there is certain probability that one of the spin 
flips lies within the region A and the other one lies within 
the region B. By projecting the evolved state by means of 
two coarse grained measurements on these regions (Middle), 
an entangled pure state between A and B can be established 
when the results of the measurements ascertain that one spin 
flip is in A and the other one in B while saying nothing about 
their exact localization within A and B (Bottom). 



where R refers to all sites in the ring. This probability 
can be interpreted as the probability of success of the 
procedure. 

Fig 3. shows the results of the simulations carried out 
in this out of equilibrium setting. We have used a ring 
with N — 24 spins and the spin flips are initially located 
at ri = 10 and r 2 = 14. As the size of measurement 
regions increase in the number of sites they contain, the 
probability of success and the Epp also increase. In the 
limit in which regions A and B equal the half of the chain, 
the probability of success Prob(\^ t ) ab) ~ 0.5 an( i the 
value of Epp slightly oscillates around 1/2. These values 
are decreasingly modulated as the size of areas are con- 
currently decreased. The time evolution induces some 
oscillatory effects in the values for entanglement, proba- 
bility and Epp for a fixed size of the measurement areas. 
For the type of system studied in this section and the 
type of projections we made, Epp turns out to be equiv- 
alent to the notion of entanglement of particles of Refs 
[HI], [IB]. In [13] this procedure has been applied to the 
continuous case of ID impenetrable bosons lying on a 
ring and authors have also shown the scaling of the Epp 
with the number of particles in the ring. An interesting 
open question is the issue of measuring the entanglement 
created by our procedure. 
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FIG. 3: Entanglement Entropy (top), probability of success, 
(middle), and Epp (bottom) of the system of two spin flips 
in the XY ring as function of time and size of the measurement 
areas A and B on the ring. Entanglement, probability of 
success and Epp scale as the size of areas increase. When the 
size of the measurement areas equal the half of the chain, Epp 
slightly oscillates around the value of 0.5. 



\$ m {p}) = — x m ({x}\{ P })\x m --- Xl ) (14) 

V m\ 

The complex valued function x is denned by 

Xm {{x}\{p}) = Yl £ ( Xfc ~ X ») det ( D {x}\{p}) 

\l<a<b<m / 

(15) 

where (-D-0}|{ P } is the m x m matrix with (j,k) entry 
Djk = exp(ixjPk)- The function x({ :z; }|{rf) is called a 
wave function and it is symmetric in x and antisymmetric 
in p. 

For some high values of the magnetic field h in (|13p the 
ground state of the ring can be expressed as the m = 2 
case of the eigenvector (fT4)) . i.e., 

1 N 

\^{p}) = ^= e(x 2 -x 1 )det(D {x}m )\{x}) (16) 

where {x} = {x\,x 2 ), \{x}) = \xxx 2 ) and det(D { X }\{ P }) 
means the determinant of the D 2x2 matrix with elements 
Di n = e^ iPnXi ' with {p} being the momenta of the two 
spin flips. We choose the same projectors Pa and Pb 
for the Pure state Entanglement Extraction as the ones 
proposed in the previous Section. The bipartite projected 
pure state l^)^ = PaPb\® 2 {p}) can be read now as, 

l*).4B = -7= e ( x k-Xj)det(D {xj!Xk}l{p} )\xj)\x k ) 

(17) 

noting that now e(xh — Xj) = +1 Vj G A, k G B. As 
in the previous section the reduced density matrix of re- 
gion A can be described by (/ \p A \j) = \ Efces A 3^ A *y k 
where Ajk = e(xu - Xj)det(D{ Xj ^ k y\{ p y). Fig 4 shows the 
values of £pp, Prob(\$>) AB ) and the entanglement of the 
projected state E(\^ab}) as the size of the measurement 
region A (B) increases. 



C. Epp in the ground state of a QMBS 

We choose to illustrate this section, the ground state of 
the ring of spins in an external field h interacting through 
the Hamiltonian 

H = H - 2hS z (13) 

where S z = E^i Noting that [H ,S Z ] = 0, the 
eigenvectors of 7i are determined by a positive integer 
m < N and a collection of real numbers {pi}™= 1 , with 
each — 7r < pi < it, called momenta [3]. For short, 
the collection {p;}™ 1 may be denoted simply by {p}. 
Thus eigenvectors of (fT3"]) are defined by the set of vectors 
|$ m M) with 



IV. GENERAL PROCEDURE FOR 
EXTRACTING PURE ENTANGLED STATES BY 
DISTANT LOCAL PROJECTIVE 
MEASUREMENTS 

In contrast to the preceeding Sections, here we elabo- 
rate on whether one can find nonzero Epp in situations 
where the properties of ground state of a QMBS (sym- 
metries, good quantum numbers) do not a priori suggest 
a heuristic choice for the suitable projectors for the pro- 
jective extraction procedure. To this end we provide a 
general prescription that allows to systematically quan- 
tify the maximally achievable Epp for a given state, when 
the distant parties A and B have access to a fixed number 
of contiguous sites. 
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which is the Schmidt decomposition of \^}. We require 
the global projectors Pj to be of direct product type 
Pj = P j A ®P B and fulfill the condition a) (Eqs. (20] ) , ([2T] ) ). 
Representation Eq. (|2"3"|) suggests the following construc- 
tion prescription 



10 12 



Y,\°i)(°i\ ( 24 ) 

The additional superscripts /i, v label different sets of in- 
dices 1^ encoding the the p, v-th. combination of rank-1 

projectors (e.g., Ic^Ko^l), leading to rank > 2 (see also 
Sec. IIII Ap . E.g. for some particular v = 3, (i = 4 and 
I 3 = {3,5}, I 4 = {2,7} the projectors in Eq. [Ml would 
become 



FIG. 4: Entanglement entropy, probability of success and £pp 
in the ground state of the ring (N = 24) with Hamiltonian 
(|13|) as the number of sites comprising the measurement re- 
gion A varies from Aa = 2 — > Aa = N/2. 



We begin with specifying, under which circumstances 
a local projective measurement on subsystems of a MBS 
can lead to a pure state. I.e. we ask for the criterion the 
projector P = PaPb must satisfy, so that 



Pp AB P^ |$X*| 



(18) 



For this purpose it is instructive to consider the reduced 
density matrix of the two separate regions A and B pab 
in its diagonal representation (Schmidt representation, 
assuming the ground state is pure) 



\l:/\®Ab\ 



(19) 



There could be two alternative ways of fulfilling 
Eq. (fl"S|) , being a) by finding a P = Pj so that 



P 3 \V) = \&) 



(20) 
(21) 



resulting in pure state |$) = |4> J ) after projection or b) 
by having a set of at least two Schmidt vectors {1^")} 
which after projection by P result in 



= |$)Va. 



(22) 



It is easy to show that one can exclude possibility b) 
by virtue of the orthogonality (^ k \^ 1 ) = 8^,1 and the 
idempotence of the projectors P 2 = P. 

In order to construct a trial set of local projectors Pa 
and Pb, we proceed by casting the vectors into the 
form 



k>i#> 



(23) 



pi 3 

A 

pi* 
B 



\4)(4\ + \4)(4\ 



We then proceed by systematically checking for the exis- 
tence of projectors piv* v ) = pw g) P B U being suitable for 
the extraction of pure states. If several such projectors 
can be found, the corresponding extracted states |$)^ A " y - ) 
will be evaluated for their entanglement content as mea- 
sured by the von Neumann entropy S^rpi^y^^ ($!))• 
£pp is then defined by the maximally achievable amount 

Pro&flS)^"))) S{Tr B {W i>lv) among these states. 
In the present notation the probabilities of obtaining the 
states read 



Prob(\<S>) 3{ ^ ] )) = \j (WIP^IW) 



(25) 



As for our procedure Eqs. (|2U)l and (|2ip express the 
only possible scenario in which pure states can be ex- 
tracted, £pp will not merely be a lower bound but quan- 
tify the optimal performance of extracting pure state en- 
tanglement by local projections. 



Application of the general procedure: 
general phase of a spin model 



£pp for a 



For the study of more general instances of Epp we now 
slightly alter the type of Hamiltonian used so far, in that 
we introduce an anisotropy causing the model to have less 
symmetries and therefore less good quantum numbers. 
The transverse XY Hamiltonian reads 

N 

HxY = -Y / ((i - 7K<+i + (i + tK4h) - 2hS * 

i=l 

(26) 

We shall explore the different ground states of TLxy for 
a particular system of N = 6 spins arranged on a ring, 
partitioned into symmetric regions A and B of two con- 
tiguous spins each that are separated by a single site. 
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FIG. 5: General procedure (top) Maximum entanglement en- 
tropy (middle) Maximum probability of extracting a pure 
state (bottom) projectively extractable pure entanglement 



For 7^0 the z-Magnetization is no longer conserved, 
leaving only the parity P z = Ylk=i a k as a commut- 
ing observable. We shall now apply our general proce- 
dure and evaluate £pp as a function of parameters h and 
7, transverse field and anisotropy respectively. Figure 
[5] shows that nonzero pure entanglement can indeed be 
extracted by local projections, even outside the previ- 
ously covered regimes, where the ground state is a su- 
perposition state of the form Eq. (|16p . In the present 
study, this corresponds to 7 = and h ~ ±0.5 . Be- 
yond these regimes our procedure reveals that pure en- 
tanglement can be extracted in a sector of the (7,^1) 
plane for which the ground states attain odd parity, i.e. 
(P z ) — — 1 , providing strong evidence that the origin of 
nonzero £pp, as detected by our general procedure, can 
be achieved by a local parity measurement. Let us ex- 



amine under which circumstances locally measuring the 
parity in regions A and B can lead to a pure entangled 
state (for simplicity restricting ourselves to the isotropic 
(7 = 0, h = 0) case, where the number of spin flips 
m = £ is conserved). An important role must be at- 
tributed to the particular partitioning we chose above, 
i.e. the complement of the measurement regions is com- 
prised of merely two sites. If the parity of regions A and 
B is even (\H), |TT))i then the complementary spins must 
be of odd parity (|TI)>|IT))) to comply with the global 
odd parity of ground state. The local measurements of 
parity must not discriminate between the particular con- 
figurations which is a necessary requirement of having 
an entangled state after projection, and is equivalent to 
the aforementioned necessity of having local projectors 
of rank higher than one. Then, measuring even parity in 

both regions A and B (Pa = P B = IT, TXT, Tl + li, IXi, il) 
is consistent with four possible spin configurations, which 
occur with equal amplitudes in the ground state, thus 
after projection the state of the whole system reads 

= ^(|U) + Ut>hfl ® (ITT)aIU) b + IU)aITT) b ) 

which was verified numerically. The state of regions A,B 
is therefore pure and entangled, which is clear by inspec- 
tion. This study therefore sheds light on the interrela- 
tion between £pp and observables that commute with 
the underlying Hamiltonian, and further have a degen- 
erate spectrum. We stress that our study has primarily 
been conducted, at the current stage, for conceptual pur- 
poses. Whether nonzero £pp exists in more realistic or 
potentially useful settings can not be said at this point. 



V. DISCUSSION 

The non equilibrium setting of Section IHI Bl shows 
qual itative similarities to the features after a local quench 
(Il9ll . [20l ]. (22I]). Calabrese and Cardy recently proposed 
[2 11 ] a simple interpretation for the observed behavior in 
terms of quasiparticle excitations emitted from the ini- 
tial state at t = that freely and semiclassicaly prop- 
agate with velocity v. If the initial state is separable (in 
the standard basis) particles that originate from different 
sites are incoherent, but pairs of particles moving to the 
left or right from a given point are highly entangled. The 
reasoning of these authors is that a point xa in region A 
will be entangled with a point xb in region B if a pair 
of entangled particles emitted at an earlier time arrive 
simultaneously at xa and xg. This picture of a causal 
cone is capable of qualitatively explaining the observed 
saturation of the block entropy of contiguous blocks of 
spins in an infinite chain after the quench [20j . It is thus 
tempting to qualitativeley interpret the results of Section 
piIBp where a saturation is observed (for a time window 
before finite size effects kick in). Our procedure of select- 
ing particular measurements outcomes assures that these 
excitations have arrived to regions A and B at the instant 
of measurement which, in agreement with the quasipar- 
ticle picture, may lead to entanglement after a time that 
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scales linearly with the distance between these regions 
due to the finite speed for the propagation of the excita- 
tions v. More importantly such a measurement does not 
discriminate particular points xa and xb within these 
regions, and therefore the quantum indeterminacy (ow- 
ing to the time evolution) is what gives us the possibility 
for establishing a pure entangled state between A and 
B. The growth and later saturation of entanglement en- 
tropy between blocks A and B can be viewed from the 
fact that due to time evolution, there is an instant t* at 
which there is a nonzero and non increasing probability 
that the excitations have arrived to all the points xa and 
xb lying within regions A and B. 

As particularly highlighted in Scc lIV Al a vital aspect of 
£pp is its connection to good quantum numbers, which 
may lead to a similar quantum indeterminacy as men- 
tioned for the dynamic case above. This can arise for 
the ground state of a QMBS if the underlying Hamilto- 
nian commutes with an observable that has a degenerate 
spectrum, like in our cases the magnetization or parity. If 
these observables can also be measured locally, and the 
local measurement outcomes are again consistent with 
more than one global state configuration (degeneracy), 
we have shown that this indeterminacy can be exploited 
in that it may lead to pure entangled states of different 
regions of the MBS after the measurement. So we are 
deliberately measuring some incomplete set of commut- 
ing observables locally, which will maintain the (quan- 
tum) uncertainty of telling in which particular state the 
subsystem resides is maintained. As these local mea- 
surements are performed on both separated regions, the 



projected state has to accomplish both local degeneracies 
simultaneously which is the origin of the nonlocal quan- 
tum correlation. Under which circumstances this quan- 
tum correlation is that of a pure state is connected rather 
subtly to an interplay of the observable being measured, 
the partitioning of the system into measurement regions 
A,B and their complement and the global properties of 
the state. 

Another interesting point to discuss is the topic related 
with the transferring of entanglement created by our pro- 
cedure to systems in which this entanglement could be 
used as a resource for quantum information tasks. In 
other words, what we would like to do is to be able 
to swap the entanglement established between regions 
A and B onto a couple of related systems 5,4 and Sb- In 
[2^ . authors propose a protocol by which the entangle- 
ment between separated regions in an ion trap could be 
transferred to a pair of ancillary ions each one of them re- 
lated with each one of the regions. Further investigations 
should clarify if the above kind of operations feasibly al- 
low to use the entanglement extracted by our procedure 
in quantum information tasks. 
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